The dark degeneracy and interacting cosmic components 
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We study some properties of the dark degeneracy, which is the fact that what we measure in 
gravitational experiments is the energy momentum tensor of the total dark sector, and any split 
into components (as in dark matter and dark energy) is arbitrary. In fact, just one dark fluid is 
necessary to obtain exactly the same cosmological and astrophysical phenomenology as the ACDM 
model. We work explicitly the first-order perturbation theory and show that beyond the linear order 
the dark degeneracy is preserved under some general assumptions. Then, we construct the dark fluid 
from a collection of interacting fluids. Finally, we try to break the degeneracy with a general class of 
couplings to baryonic matter. Nonetheless, we show that these interactions can also be understood 
in the context of the ACDM model as between dark matter and baryons. For this last investigation 
we choose two independent parametrizations for the interactions, one inspired by electromagnetism 
and the other by chameleon theories. Then, we constrain them with a joint analysis of CMB and 
supernovae observational data. 

PACS numbers: 98.80.-k, 98.80.Es, 95.35. +d,95. 36. +x 
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I. INTRODUCTION 

Current cosmological observations indicate that about 
96% of the total energy content of our Universe is made 
of yet unknown dark components and only 4% is made 
of particles of the standard model. Among them, there 
are precision measurements of anisotropies in the Cosmic 
Microwave Background (CMB) radiation [3-Q, baryon 
acoustic oscillations [3, [f| , and Type la Supernovae p-ll] . 
For a recent review of the nowadays status of cosmology 
see, e.g., @. 

Usually this dark fluid is separated into two compo- 
nents: a clustering dark matter piece responsible for 
forming structure in the Universe and a nonclustering 
dark energy with a negative pressure that is responsible 
for the current accelerated cosmic expansion. 

Moreover, in the concordance model of cosmology -the 
so-called Lambda Cold Dark Matter (ACDM) model- the 
dark energy is of geometric nature and enters as a con- 
stant at the Lagrangian level of Einstein's gravitational 
theory. For gravity matters, the cosmological constant is 
indistinguishable from the vacuum contribution of quan- 
tum fields, and here appears an intriguing problem, the 
sum of both contributions is about 1 over 10 120 times 
the latter 10]. A second conceptual problem with the 
ACDM model comes out when one considers the ratio of 
dark energy over dark matter energy densities, which in 
spite of growing with the third power of the scale fac- 
tor of the Universe, its value today is of order one. This 
is the so-called coincidence problem. As a consequence, 
a plethora of alternative proposals has appeared in the 
literature; as a short sample see (TT1 - [20| [ . 



In this paper we argue that the separation of the dark 
fluid into dark matter and dark energy is arbitrary and 
is only favored for historical reasons and computational 
simplicity. After all, we define the dark fluid as our lack 
of knowledge as 1 
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where comes from the observed geometry of the Uni- 
verse and T°^ s from its observed energy content. In fact, 
the dark sector could be composed by a large zoo of par- 
ticles with complicated interactions between them. Or, it 
could be even just one exotic unknown dark fluid. This 
property has been called dark degeneracy by M. Kunz 
[HI; see also (HHH)]. 

In this paper we assume an astrophysical perspective 
to the dark fluid by defining it as a barotropic fluid with 
speed of sound equal to zero, as in [29l l30ll and more 
recently in [31(; for similar approaches see [33 - 1341 ] . 

Making the speed of sound equal to zero ensures that 
dark fluid energy density perturbations will grow at all 
scales, but at the same time, it allows the fluid to have a 
nonzero pressure. This is quite in contrast with canonical 
scalar fields as quintessence, for which the speed of its 
perturbations is equal to one, a feature that makes dark 
energy perturbations to be quickly damped, so they do 
not grow inside the horizon. 

Interactions within the dark sector has been studied 
largely in the literature [35l - f43| . mainly as a mechanism 
to solve the coincidence problem. On the other hand, in- 
teractions between dark matter and the standard model 
of particles are expected, and have been studied also on 
several occasions. In fact, the weakly interacting massive 
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We use natural units, h = 1 and c = 1, otherwise is stated. 
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particles (WIMPs) paradigm has emerged as the predom- 
inant scenario to solve the missing mass problem [44r - 
HU. Nevertheless, other alternatives have been proposed. 
Among them, special attention has been attracted by the 
strong interacting dark matter scenarios |49l - l54j , in which 
dark matter interacts with itself and with bary ons [Hoi ] 
through the strong force. It has been shown [49ll55j that 
this alternative can alleviate the cuspy halos and over- 
density of substructure problems that appear in N-body 
simulations based on the standard ACDM model [561 ] . 

In this paper we develop a general class of couplings 
of the dark fluid to the standard model of particles, and 
show that they can be also understood as interactions of 
dark matter to baryons. Then, we impose constraints by 
using CMB anisotropics and supernovae observations. 

This work is organized as follows: In Sec. II we de- 
fine the dark fluid from the properties of the dark matter 
itself. In Sec. Ill we work out the cosmological back- 
ground solutions for the dark fluid and show that they 
are identical to the ACDM model ones, leading to the 
dark degeneracy. In Sec. IV we show that the dark de- 
generacy is preserved under some assumptions when one 
goes beyond zero order in cosmological perturbation the- 
ory. In Sec. V we work a multifluid description of the 
dark sector, allowing interaction among its components. 
In Sec. VI we extend these interactions to baryons to try 
to break the dark degeneracy. Finally, in Sec. VII, we 
present our conclusions. 



II. THE SOUND SPEED OF THE DARK FLUID 

Considering Newtonian gravity for the moment, the 
ovcrdcnsitics of matter in an expanding Universe follow 
the equation in Fourier space 



5" + 2HS' 
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AirGp) 6 = 0, 
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where /c p h ys is a physical wavelength (contrary to the co- 
moving wavelength which we will use in the following 
sections) related to a physical length scale by /c p hy S = 
27r/iph ys , c 2 is the speed of sound of the fluid under con- 
sideration and prime means derivative with respect to 
cosmic time. In astrophysical scales it is safe to set H = 
and then, the scale factor is equal to a constant that we 
make equal to one for this discussion matters. It follows 
that there is a threshold scale, the Jeans length, given by 



(3) 



for which perturbations with physical length above it, 
^phys > I J, grow by gravitational collapse, and perturba- 
tions with ^phys < I j develop acoustic oscillations. 

We expect to have dark matter structure at a wide 
range of scales: from the largest cosmological structure 
to galaxies. In fact even dwarf galaxies need to have 
dark matter halos in order to stabilize their disks and to 



account for the necessary gravity sources to obtain the 
observed flat rotation curves which follow the luminous 
matter in there. The only way to guarantee that dark 
matter perturbations grow at all scales is demanding that 
the Jeans length be equal to zero. 

For barotropic fluids the adiabatic speed of sound co- 
incides with the speed at which perturbations propagate 
in the fluid. We define the dark fluid as a barotropic fluid 
with an adiabatic speed of sound equal to zero, 



0. 



(4) 



and that, at a first approximation, does not interact with 
particles of the standard model. 2 Without lost of gener- 
ality we can write its equation of state as 



P(p) = w(p)p 



(5) 



where w is its equation of state parameter and is a 
function of the energy density p only. Thus we have 
c 2 = (dP/dp) s = dP/dp. Using Eq. (g]) we obtain that 
the equation of state parameter is solved by the equation 
pdw/dp + w = 0, whose solution is 



C 
P 



(6) 



with C a constant and the minus sign is set for later 
convenience. This means that the pressure is a constant 



P = Po 



-C. 



(7) 



Astrophysical observations constrain this value to be very 
small, |P| <C pa, where pa is the energy density of typical 
astrophysical scales where dark matter has been detected. 
Usually, it is assumed that dark matter is pressureless, 
but this is by not means necessary, for instance it could 
be the case that |P| ~ p c0l where p c0 is a typical cosmo- 
logical energy density scale at present, without getting 
in contradiction with observations. In fact, this is the 
entrance that leads us to consider the dark fluid to be 
dark energy as well as dark matter. 



III. DARK FLUID AS DARK ENERGY 

Now, let us consider a Friedmann-Robertson- Walker 
description of the Universe at very large scales, filled 
with standard model particles (6, 7, ...) and with the 
above-defined dark fluid, from now on labeled by d. The 
evolution equations of such a Universe are 



2 The barotropic condition dismisses, for instance, the possibility 
of scalar fields for which although the speed of propagation of its 
perturbations is equal to the speed of light, this does not coincide 
with its adiabatic speed of sound, which indeed could be zero. 
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H 2 = 



8ttG 



{Pd + Pb + Pi), 



(8) 



behavior of the dark fluid as dark matter in astrophysical 
scenarios. 

The Friedmann equation ([8j becomes 



and 



P'd 



p' b + 3H Pb = 0, 
p' 7 + AH Pl = 0, 

f 3^(1 + w d )p d = 0, 
factor. 



(9) 



(10) 



(11) 



Hubble 
~ 3 and 



Equations 
' 4 , where 



where H = a' /a is the 

and (|10l) give pt, = Pbo^" and P-y = P-yO a ", wtiere a 
subindex means that the quantity under consideration 
is evaluated at present time, and we have normalized the 
scale factor to be equal to one today, do = 1- Integration 
of Eq. (fTTj) gives, using Eq. (j6|), 



H 2 = 



IttG 



fCp 



PdO 



1 + K 1 + K 



d0 a- 3 



Pboa 



-3 



p l0 a 



(16) 

This is the same evolution equation for the scale factor as 
in the ACDM model. This is not an accident, if one as- 
sumes ACDM as the valid model and considers the total 
energy density of the dark components, pr — Pdm + Pa, 
the total equation of state parameter, wt, defined by 



Wt 



_ Eq W aPa 



(17) 



where the subindex a runs over dark matter (DM) and 
cosmological constant (A), is given by 



Pd 



PdO 

1 + K 



(12) 



where we have defined the constant K, = (pdo — C)/C. 
This expression is what we expect for the evolution of a 
unified fluid: a piece that redshifts with the scale factor 
as a~ 3 , as a dark matter component does, and a piece 
that remains constant, as vacuum energy. This result has 
been found elsewhere in the literature, as in [34j . which 
investigated the ACDM limit of a Chapliying gas, or as 
in the study of barotropic fluids with constant speed of 
sound |29j - [33 |. In fact, using the language of E. Linder 
and R. Scherrer [32j], in our case the term proportional 
to a -3 is the aether piece of a barotropic fluid. 

Now, the equation of state parameter of the dark fluid 
fl6l) becomes 



1 



1 + /Ca- 3 ■ 



(13) 



We want to stress that the energy density of the dark 
fluid is proportional to the inverse of its equation of state 
parameter 



Pd = 



PdO 



(1 + K) w d ' 



(14) 



and that its pressure, expressed in terms of the constant 
K instead of C, is 



P,, 



PdO 

1 + K.' 



(15) 



In order to ensure the positivity of the energy density at 
all times, K must be a positive number. This implies that 
the pressure is negative, a quality that allows the dark 
fluid to accelerate the Universe, and as we have outlined 
in the last section it could take values of the order of 
the critical density (~ 3Hq /8itG) without affecting the 



wt 



1 + %M a -3' 



(18) 



where 17dm = 87rGpDMo/3i?o an d ^ao = A/3Hq. 
(Clearly, the equation pr = — 3"H(1 + wt)pt is accom- 
plished for the total energy density.) Then, comparing 
these results to Eqs. (|13l) and (fTo]) . we note that under 
the identifications 



and 



K = 



9. 



DM 



DM 



(19) 



(20) 



where fl d — 8irGpdo/3HQ, the resulting evolution cos- 
mology in both models are exactly the same, at least at 
the background level. 

This property has been called dark degeneracy by M. 
Kunz in [21|. In fact, it is more general than for the sin- 
gle fluid case worked here: any collection of fluids whose 
total equation of state parameter is equal to Eq. (fl8| and 
that do not interact with baryons and photons will be- 
have exactly as the composed dark matter-cosmological 
constant fluid, leading to a degeneracy with the ACDM 
model. 



IV. BEYOND THE BACKGROUND 

In this section we explicitly show that the degener- 
acy is preserved when one goes beyond the homogeneous 
and isotropic cosmology, but only under some general as- 
sumptions that have been taken for granted in previous 
works. We divide the discussion into linear and higher 
orders in cosmological perturbation theory. 
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A. Linear order 

Let us consider cosmological perturbation theory in the 
conformal Newtonian gauge, the metric is given by (for 
details see |57j ) 

ds 2 = a 2 (r) [ - (1 + 2*)dr 2 + (1 - 2$)<%cfoW ] , (21) 

where t is the conformal time related to the cosmic time 
by dt — adr. As usual, we define the matter pertur- 
bation variables through the expressions of the energy 
momentum tensor 

T° = -p{l + S), (22) 
T\ = -{p + P)v\ (23) 
T 3 = P((l + ^)^+n t J ), (24) 

where Wj is the anisotropic (traceless) stress tensor. The 
energy density p and the pressure P denotes background 
quantities and are functions of the conformal time only. 
The vector v l is called the peculiar velocity and is related 
to the four velocity w M of the fluid by the relation v % = 
u l /u°. In Fourier space we define the velocity 6 — ikiV 1 
and the scalar anisotropic stress a — 2kikjIP : 'w/3(l + 
w). We use the flat space metric (Sij) to raise and lower 
indices of intrinsic space geometrical objects like v l and 

The hydrodynamical equations for a general fluid are 

m 



k 2 <P = -AirGa 2 ^ p^, (31) 

i 

and 

le 1 (9-9) = 12irGa 2 ^2(p i +P i )a i (32) 

i 

where the sum runs over all fluid contributions and 

A i = 6 i + 3H(l + w i )-± (33) 

is the rest fluid energy density (58j . 

To solve these equations, we need to add information 
about the nature of the dark fluid. The barotropic con- 
dition implies that SP = c 2 5p and thus SPd = 0, and 
because it is a perfect fluid, the anisotropic stress van- 
ishes, HJj = 0. Accordingly, the space-space components 
of the perturbed energy momentum tensor are equal to 
zero, 8T¥ = 0. Therefore, in the right-hand side (rhs) 
of Eq. (|2T|) the last term vanishes, and in the rhs of Eq. 
(|28f only the two first terms survive. If there are only 
baryons and dark fluid the two gravitational potentials 
are equal, $ = 

Figure 1 shows the evolution of a mode k — 0.05 Mpc -1 
of the perturbations variables. The dashed lines shows 
results for the ACDM model, for which we have to use 
the dark matter perturbations equations, 



A = -d + "■)((-> - :-!<I>) - :-YH [ ^ - w ) ,v j2.",i 



-H(l-3w)d 
-k 2 ^ - k 2 a, 



SP/Sp, 2 



1 + w 1 + w 



k 2 S 



(26) 



where 5P = Pttl, Sp = p5, % = a/a and a dot means 
derivative with respect to conformal time. For the dark 
fluid case, c 2 = Pd/f>d — Wd — Wd/3H(1 + u>d) = and 
these equations become 



S d = -(l + w d )(6d-m + 3Hw d 6d 

-m & -f±8 dl (27) 

OPd 

9 d = -H6d + k 2 ^+ 6 -p^k 2 Sd-k 2 a d , (28) 

1 + Wd 

while for baryons after recombination (when the coupling 
to photons can be safety neglected) 



S b = -0 fc + 3$, 
6 b = -H9 b + k 2 V. 



(29) 
(30) 



The fluid equations are supplemented with the Einstein's 
equations 



0dm = -ne DM + k 2i & , 



(34) 
(35) 



instead of Eqs. ([27]) and (|28|). 

In evolving the perturbations for both models, ACDM 
and dark fluid, we have imposed on the initial conditions 
the relations 5d(n) = apT>M{Ti)&T>M{n) / Pd(n) for the 
density contrasts, and 9d(n) — ap D M(Tj)^DM('r l )/(l + 
Wd)pd(Ti) for the velocities, and we let a to take different 
values. These initial conditions are given at an initial 
time well after recombination, so we can use Eqs. (|2"9")l 
and (|30f and the relation i5dm — $b holds. 

We note that the evolution of the density contrast of 
baryons for the specific mode we have chosen is indistin- 
guishable for both models if we take a = 1. Then, al- 
though the cosmological observable is the baryonic mat- 
ter power spectrum which includes a wide range of wave- 
lengths, Fig. 1 suggests that indeed the results are the 
same for both models for any wavelength, as we show 
below. 

In the cosmological context, imposing these two initial 
conditions is equivalent to demand that at first order in 
perturbation theory, the time-time and time-space com- 
ponents of the perturbed energy momentum tensor of 
the dark fluid and ACDM models are equal at the given 
initial time r^. It is straightforward to show that both 
conditions, in the case a = 1, will be preserved at all 
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FIG. 1: Evolution of perturbation variables for a mode 
k = 0.05 Mpc -1 . Solid curves are obtained from the dark 
fluid model for different values of the parameter a in the ini- 
tial conditions 5d(Tj) = apY>M(Ti)5oM(Ti) / pd{Ti) . a takes val- 
ues from 0.8 to 1.2. Dashed curves are for the ACDM model 
variables. The panels show: (a) Gravitational potential 
(b) Baryonic density contrast 5b- (c) Dark fluid (solid lines) 
and dark matter (dashed line) density contrasts, Sd and <5dm- 
(d) Dark fluid (solid lines) and dark matter (dashed line) ve- 
locities, 8d and Bum- The solutions for the case a = 1 are 
depicted with the thick (red) lines, which for panels (a), (b), 
and (d) coincide with the dashed lines. 



times. This means that in obtaining the a 
in Fig. 1 we have imposed that 



and 



PdSd = Pdm<$dm, 



Pd(l + Wd)0d = Pdm#dm- 



1 solutions 



(36) 



(37) 



From now on, we will demand Eqs. (|36|) and (|37|) to hold 
for all considered fluids, staying in degeneracy with the 
ACDM model at first order in cosmological perturbation 
theory. Later on, in Sec. VI, we will try to break this 
degeneracy, but through couplings of the dark fluid to 
baryonic matter. 

From the results of the previous section it is straight- 
forward to see that pDM/pd = 1 + Wd, then Eq. (f3"7| 
implies that 0dm = Od, which explains the behavior of 
the curves in Fig. Id. Moreover, from Eq. (f3"Tj) it follows 
that the gravitational potentials $ are the same for both 
models, so the behavior in Fig. la. 

Finally, inserting Eqs. ([31) and into Eqs. ([51)1 
and (|35|) . respectively, and using Eq. (fl4|) . we obtain 
that the hydrodynamical perturbation equations for the 
dark fluid are 



5 d = -(1 + Wd)(0 d 

o d = -ue d + kH. 



3$) + ZUwd&d, 



(38) 
(39) 



These equations are equal to (|27|) and (|28|) for a 
barotropic perfect fluid, confirming the dark degeneracy 
at first order level. 

It is interesting to note the behavior of the density con- 
trast of the dark fluid (Fig. lc), which initially grows as a 
dark matter component to later decay asymptotically to 
zero. In fact, this behavior is natural for all wavelengths 
modes as can be seen from Eqs. (|27|l and ([280: At late 
times, when Wd tends to -1, the &d source becomes inde- 
pendent of Od and is a negative quantity which tends to 
zero. Meanwhile at early times, for w d — > 0, Hwd — > 0, 
and the equations become equal to the cold dark matter 
evolution equations, ([34]) and ([35]) . 



B. Beyond linear order 

To go beyond the linear order, let us make pertur- 
bation expansions to the dark fluid and ACDM energy 
momentum tensors about the (zero order) background 
cosmological fluids as 



T — IT(0) i 7-l(l) , rp(2) 1 _ § 



(40) 



If the total energy momentum tensors of both models 
are equal (T£ v — T^ CDM ), clearly each of the terms in 



the expansion will be equal as well {T^^ = T) 



- n ACDM(i)N | 
flu j 
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This argument is outlined in [2l| to argue that the degen- 
eracy is preserved at all orders. Certainly, this is correct. 

Nonetheless, we want to stress a different approach: 
We are affected gravitationally by the total energy mo- 
mentum tensor, but usually when comparing observa- 
tions to models we expand it as in Eq. (j40|) . and after 
this we assign values to each of the pieces. The fact that 
both energy momentum tensors are equal, say at zero 
order, does not imply that they will be equal at first or- 
der. In this situation, equations such as (|3"o) and (|57|) 
are conditions of the theory and not consequences of it, 
and if not imposed, the degeneracy is broken, as seen in 
Fig.l for a^l. 

This approach, which we will adopt, is very close to 
that followed in some papers which investigate modified 
gravity theories that are indistinguishable in the back- 
ground from the ACDM model but could differ in linear 
order, leaving imprints that are parametrized to compare 
to current and future observations (see [H^] and references 
therein). In fact, if some of these imprints were detected, 
they could be a consequence of interactions within the 
dark sector (or even with baryons, as those we study in 
Sec. VI), instead of a deviation of General Relativity. 

Finally, there is a subtle point that is pertinent to as- 
sert at this moment. The fact that the Universe is so 
smooth at very large scales is what allows us to expand 
the energy momentum tensor of the cosmic fluids as in 
Eq. (|40|) . The zero-order term is in fact a spatial average 

in hypersurfaces of constant cosmic time, T$) = 
and thereafter we assume that general relativity holds 
in the form G flv (g^) = (T^ u ). This is not true in gen- 
eral due to the nonlinear character of gravity; it is well 
known that backreaction effects contributes to this last 
equation. For recent work in the averaging procedure in 
cosmology and its caveats see (60j . 

For the specific case that concerns us there is a re- 
lated issue. Let us consider a fluid with equation of state 
P = P{p) that leads to a unified description of the dark 
sector. It has been noted that it is not necessarily true 
that (P) = P((p)) [61]; in fact higher-order corrections 
enter into this equation. Therefore, when some scale 
grows and leaves the linear regime, the naive averaging 
procedure on the equation of state is no longer valid. This 
nonlinear effect has been investigated in [23 . l6ll |62| , and 
it was shown that instabilities are expected to occur in 
unified dark models, even on large cosmological scales. 
Two exceptions are the cases in which the involved equa- 
tions of state are P = constant and P = wp, with w 
a constant. The former is the case of the dark fluid, 
and then, it does not suffer for this averaging problem. 
Nonetheless, it is worth emphasizing that in models that 
depart -although indistinguishable by current observa- 
tions when calculations are performed using the average 
procedure- from the dark fluid or from the ACDM mod- 
els these effects must be considered. 



V. MULTIPLE INTERACTING DARK FIELDS 

From the last section it is clear that the dark fluid, al- 
though not necessarily, could be the sum of a dark matter 
and a cosmological constant components. Nonetheless, if 
the interaction to the particles of the standard model 
is only gravitational (which is the ultimate definition of 
dark), the nature of the dark fluid is fundamentally im- 
possible to elucidate, because of the universality of this 
force. 

In this section we explore the possibility that the dark 
fluid is composed of a collection of dark components with 
possible interactions between them. The equations that 
we obtain will be generalized in the next section to in- 
clude interactions to baryons. Let us consider the Ein- 
stein field equations, 



G^ = 87rG\^2Tr + TZj, (41) 

where the subindex a labels the different dark compo- 
nents. If we forbid dark components to interact with 
standard model particles, the Bianchi identities imply 

^m^sm = and 

VplT = Q v a , (42) 

where the energy momentum transfer vectors, Q^, obey 
the constraint 

E^=0. (43) 

a 

Considering for the moment the homogeneous and 
isotropic cosmology, the continuity equation for each fluid 
is 

pa + 3H{l + w a ) Pa = q a , (44) 

where q a denote the energy transfer between different 
dark components. This implies that at zero order, (ig- 
noring the subindex a's for the moment) 

Q"C°> = l(g,0), (45) 

if we define Q = yf—g^vQ^Q", then Q(°> = q/a and 
qu(o) = a -i(Q(0) 5 o). Now, following [§3, [H, we split 
the interaction vector as 

Q v = Quu + U, (46) 

where the momentum transfer, f v , is first order and is 
orthogonal to the four velocity j v u v — 0. Thus, under 
spatial rotations Q transforms as a scalar and /* as a 
vector. Accordingly, we decompose f l as 

fi = f,i + *i (47) 
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with / a scalar and e, a transverse vector, e i i = 0, and 
we define the perturbation <5Q through 



and 



1 



Q = QW+6Q=-(q + 5q), 
a 



(48) 



where in the last equality we defined 5q = aSQ. From 
(l46l) it follows that up to first order 



Q° = _(g(l + 

/l^ /T 1 



(49) 
(50) 



Finally, from Eq. (|43|) one obtains (after restoring a's) 
the constriction equations 



5>« = o, 



(51) 



and a 2 £ Q* = E„W + /a + 4). which b Y takin g 
the divergence and going to Fourier space become 



Ma + k 2 f a ) = 0. 



(52) 



Now, to first order in perturbation theory, the divergence 
of the energy momentum tensor -for each fluid, omitting 
the subindex a- is 



V M T"° = 1(1 + 6 -2*){p + 3U(l+w)p) 



5 + (l + to)(9y-3$) 



:;?/ [ ^ - to I S 



(53) 



= — djip j + — { P + m(i + w)p)i 



p + p 



pw 



H(l+w)p 



1 + w 



1 + w 



(54) 



Note that if the interaction vanishes, both equations re- 
duce to the usual ones. We take the divergence of Eq. 
([54")) to isolate the scalar mode perturbations, then Eq. 
(|42|) implies 



5 a + (I + w a )(9 a - 3$) + 
3n( 5 -f^-w a )s a + ?±(S a -*)- S -*i =0(55) 

\OPa J Pa Pa 



9 a + U 1 - 3u; Q 



q a w a 



n(l+Wa) P a ° 



+ - 



SPa/Spa 



l+W a 1 + Wa 



p a (l + U) a ) 



0. 



(56) 



where we have restored the index a's and back to Fourier 
space. Also we have used the decomposition for IP- 7 in 
tensor, vector and scalar pieces. Then when hit with 
kikj only the scalar anisotropic stress survives and the 
first term of the rhs of Eq. ^ becomes 2fc 2 n>'P/3 = 
(p + P)k 2 a. 

Defining the total energy density of these dark fluids 
as pt — ^2 p a , it follows that the total density contrast, 
5t, is the weighted sum of the individual components 



Pt 



and the total velocity is given by 



(57) 



7 T = 



rX>„(l + U>«)0«, (58) 

Pt(1 + wt) ^ 



where wt is given by Eq. (|18[) . and the index a runs 
over all dark components. Using these two last equations 
it is straightforward to show that the hydrodynamical 
equations for the total fluid are 



S T + (l + w T )(0T~'S^) + m 



SP T 



\§1>i 



WT )5 T = (59) 



and 



6 T + H(l-3w T )0 T + 

5P T /5p T ,2, 
— k 0t 

1 + WT 



wt 



1 + wt 



T 



k 2 ^ 



k 2 <JT 



0, (60) 



where 8Pt = ^26P a - These are exactly the equa- 
tions of a single noninteracting fluid, Eqs. (1251) and 
(f2l)|). The adiabatic speed of sound of the total fluid is 
c 2 T = Pt/pt = J2PaC 2 sa /PT- Thus, by forcing the con- 
dition c 2 T = we obtain that wt = Wd, and if each dark 
component has positive energy (p a > 0), it also implies 
that c 2 a = and SP a = 0. Therefore, we obtain that 
the hydrodynamical equations for the total fluid pertur- 
bations become the same as the dark fluid perturbations 
(Eqs. (jMJ) and 433)), under the substitution T -> d. 
Thus, the composed fluid we have constructed is just the 
dark fluid. 

If one considers the unphysical case in which not all 
dark components have positive energy density, the speed 



8 



of sound of some components can be different from zero, 
and the total pressure perturbation is equal to 



5Pt = (? b t&Pt ~ 



1 



6Hp T 

a,b 



where 



Sab 



S a 



1 



Paf>b(c 2 sa ~ C 2 sb )S ab , (61) 



(62) 



1 



Wb 



is the relative entropy perturbation between fluids a and 
b 65[. Therefore, to obtain the dark fluid in the case in 
which some components have negative energy density we 
have to impose the pressure perturbations given by Eq. 
(f6~Tj) to be equal to zero. This condition also implies that 
the total fluid is barotropic. 



VI. INTERACTIONS TO THE STANDARD 
MODEL OF PARTICLES. 
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FIG. 2: The CMB power spectrum considering different val- 
ues of the interaction parameter Sj in units of 10~ 6 ot. The 
other parameters are fixed. 



In this section we study dark fluid couplings to bary- 
onic matter. We consider models in which the back- 
ground cosmology is the same as in the ACDM model, 
accordingly we do not allow energy transfer (q = 0) be- 
tween the cosmic components. Nevertheless, we allow a 
momentum transfer different from zero. Thus, the inter- 
actions affect the fluids only at first order in perturbation 
theory. 

In this work we will not consider interactions to elec- 
tromagnetism. This is not only for simplicity, many 
theoretical models present conformal couplings, as the 
chameleon theories 
sor gravity theories 



6000 



66, 67[ (and in general, scalar ten- 
68|), or even direct co upl i ngs to the 
trace of the energy momentum tensor p7l l69l - [72| . Also, 
this is expected in scenarios like the strong interacting 
dark matter (4t| where the couplings are through 
the strong force, for which photons are chargeless. 

In the previous section we found the hydrodynamical 
equations for a system of coupled dark components. It is 
straightforward to generalize those results to a coupling 
between the dark fluid and baryonic matter. The conser- 
vation equations for the perturbations are, for the dark 
fluid, 

S d = -(l + ™ d )(0 d -3$)+37W d + — (63) 

Pd 

= -ne d + k 2 ^ — fd - , (64) 



p d {l +w d )' 



and for baryons 



Pb 



9 b = -H9 b + k z V + c 2 sb k 2 6 b 



k 2 f b 



p b (l + w b ) 



(65) 
(66) 



These equations are supplemented by the constrictions 
(|5ip and (fS"2"]) . The interactions to electromagnetism are 



Cj 3000 



1000 




FIG. 3: The CMB power spectrum considering different val- 
ues of the interaction parameter £// in units of ctt- The other 
parameters are fixed. 



considered, but not shown in Eq. (IrJBl . Note that we have 
not considered a term c 2 d k 2 Sd in Eq. (|64[) . in agreement 
with the definition of dark fluid. 

In the absence of a fundamental theory we parametrize 
the coupling. But before do it, let us make a compari- 
son to electromagnetism. Under the formalism we devel- 
oped in the last section, the hydrodynamical perturba- 
tion equations for Thomson scattering (neglecting photon 
moments beyond the quadrupole) can be obtained if we 
choose <5g 7 = and k 2 f 1 = — p 7 (l + Wj)ax e n e aTc(6 — 
6* 7 ), where <tt = 6.65 x 10 _25 cm 2 is the Thomson scatter- 
ing cross section, n e is the number density of electrons, 
x e is the ionization fraction, and c = 1 is the velocity of 
light 

Using the analogy to electromagnetism we choose the 
interactions terms to be 
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FIG. 4: (color online) Marginalized probability for the complete set of parameters. Solid lines (red) are for Model A, dotted 
lines (blue) for Model B, and double-dotted (black) lines for Model C. The data used are the WMAP seven-year results, Union 
2 data set supernovae compilation and a prior of HST on the Hubble constant. 



and 



6g d = 0, (67) 



f d = Pd {l + w d )an d Hi{e b -e d )/k 2 , (68) 



where the parameter Ej has units of area times veloc- 
ity, or thermalized cross section (<jv), which we identify 
with some, unknown, fundamental interaction. n d is the 
number density of dark particles that we set equal to 



n d 



PdO 

m p a 3 



(69) 



where we use the proton mass, m p = 0.938 GeV, as an 
arbitrary mass scale and p d o is the energy density of the 
dark fluid evaluated today. Here, we have not an analo- 
gous to the ionization fraction, in empathy to universal 
interactions. 

The constrictions fa = —fa and 5q d — —Sqt, = 0, 
enter into the baryons equations. Because of the re- 
lation pdm = (1 + Wdjpd t ne interaction term goes as 



(dd — 8b)a~ 2 ■ Then, the greatest deviation from the stan- 
dard model of cosmology will come from the early Uni- 
verse. Accordingly we expect the interaction to be tightly 
bounded from CMB anisotropics tests. We have modi- 
fied the public available code CAMB 73] to account for 
the interaction. 

In Fig. 2 we show the angular power spectrum for the 
CMB radiation. There several curves are shown for differ- 
ent values of Ej, which is expressed in units of 10 -6 times 
the Thomson scattering cross section times the speed 
of light (c = 1). We note that the greatest differences 
with the ACDM model comes at high multipoles, which 
is not surprising because the interaction quickly decays 
and only modes which have entered into the horizon at 
early times in the Universe were affected by it. 

On the other hand, it is well known that interactions 
between the dark sector and baryons are tightly con- 
strained from equivalence principle and solar system tests 
74j . Nonetheless, such interactions could be originated 
by general predictions of string theory [7^, Ij^l . leading to 
the necessity of some screening mechanism [66ll67ll77ll78j] 
to evade these experimental constraints. In chameleon 
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TABLE I: Summary of constraints. The upper panel contains the parameter spaces explored with MCMC for each one of 
the three models. The bottom panel contains derived parameters. The data used are the WMAP seven-year data, Union 2 
compilation and HST. 
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theories [(3(| |67j], the range of this interaction depends 
on the energy density of the surrounding medium, lead- 
ing to long-range forces when the density is low, and to 
short-range forces when it is high. Inspired by these theo- 
ries, we give a k dependence to the hydrodynamical equa- 
tions, such that the interactions respond to an effective 
wavenumber given by k 2 s = k 2 g(p), where g is a mono- 
tonic growing function of the ambient energy density. By 
noting that p falls with the scale factor, we choose a po- 
tential law for the wavelength number, k 2 s = k 2 /a n . In 
the following we specialize to the case n = 2. We use the 
substitution k — >• k c g in Eq. and obtain 3 



the effective interaction remains constant while the scale 
factor grows and all modes feel it when they enter the 
horizon. 

It is possible to treat both parametrizations together 
if we choose 

Jd = + w d ) = (9 b - B d )/k . (71) 

The hydrodynamical equations become 



f d = Pd (i + w d )z n ^(e b -d d )/k 2 . (70) 

m p 

Clearly, we can also understand this interaction as a de- 
pendence of the cross section on the ambient energy den- 
sity. In Fig. 3 we show the CMB angular power spec- 
trum for different values of the parameter £//. We note 
that, in remarkable contrast to the £/ interaction, here 
all scales are affected in a similar way. This is because 



3 Chameleon theories are modifications of gravity and their per- 
turbation equations are different from those used here; it is out 
of the scope of this work to treat the precise equations of the 
chameleons. For such a treatment see l79t l8dH . 



5 d = -(l + w d ){9 d -3<l>) + mw d 5 d , 

e d = -ne. + kH- ^ + ^f^ 



and for baryons 

5 b = -0 6 + 3$, 

9 b = -H6 b + kH + c 2 sb k 2 5 b 



(72) 

(9 b - 9 d ),{73) 



(74) 



Pd 
Pb 



(l + w d ) 



(Sj ± T, H a 2 )p d0 



(75) 



If one assumes the ACDM decomposition of the dark fluid 
in dark matter and cosmological constant, it is easy to 
see that instead of Eqs. (|63l) and (|64l) . the following 
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FIG. 7: Contour confidence intervals for ilth 2 vs J7dm/i 2 at 
68% and 95% c.l. Solid (red) lines are for Model A, dot-dashed 
(blue) for Model B, and thick dashed (black) for Model C. 



FIG. 5: Contour confidence intervals for Ej vs Ej/ at 68% 
and 95% c.l. The shading shows the mean likelihood of the 
samples 





FIG. 6: 2D Posterior confidence intervals for fij/i 2 vs E/ 
at 68% and 95% c.l. Left panel: Considering both interac- 
tions (model A). Right panel: considering only interaction 
E/ (model B). The shading shows the mean likelihood of the 
samples 
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FIG. 8: Contour confidence intervals for n s vs A s at 68% and 
95% c.l. Solid (red) lines are for Model A, dot-dashed (blue) 
for Model B, and thick dashed (black) for Model C. 



equations govern the evolution of the dark matter per- 
turbation 



Sum = -0dm + 3$ + ^1, (76) 

PDM 

0bm = -W DM + fc 2 f+^, (77) 

PDM 

where we have used the conditions (|56")l and (I57|) . The 
transfer energy and momentum terms are related by 



Sqdm = 8q d and /dm = fd- (78) 

Thus, although the degeneracy with the ACDM has been 
broken at first order in perturbation theory, there exist 
degeneracies to other models, as in this case to ACDM 
with the same interactions, which means that the general 
class of interactions given by Eqs. and does not 
help us to elucidate the actual decomposition of the dark 
fluid (if it exists). Accordingly we can treat the above 
described couplings as interactions between dark matter 
and baryonic matter, without loss of generality, as we 
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FIG. 9: CMB power spectrum for the mean values estimated 
with the MCMC analysis by CosmoMC for the three models. 
The reported difference is A = 1(1 + 1) (cf° dcl c - Q). 

will do for numerical purposes 4 . Thereafter, if desired, 
we can go back and forth between both models' results 
using Eqs. ([19]) and (|20|) . 

To constrain the interactions, we perform a 
Monte Carlo Markov Chain (MCMC) analy- 
sis over the nine-parameter space (Model A) 
{Of,/i 2 , fWfr 2 , 6, T,n M , log A s , A sz , using the 

code CosmoMC [81| . 9 is the ratio of the sound horizon 
to the angular diameter distance at recombination, r is 
the reionization optical depth, n s is the spectral index 
of the primordial scalar perturbations and A s is its 
amplitude at a pivot scale of ko = 0.05 Mpc -1 . 

We have imposed flat priors on the two interaction 
parameters: < Ej < 10~ 7 x ctt and —11 x <jt < 
Yin < 10 x (Tt- For the CMB anisotropies and polariza- 
tion data we used the Wilkinson Microwave Anisotropy 
Probe (WMAP) seven-year observations results [3|. For 
the joint analysis we use also Hubble Space Telescope 
measurements (HST) [82l | to impose a Gaussian prior on 
the Hubble constant today of Ho = 74 ± 3.6km/s/Mpc, 
and the supernovae type la Union 2 data set compilation 
by the Supernovae Cosmology Project 0. 

We also study two other models: Model B, only 
considering the interaction Ej, it has an eight-parameter 
space {fifc/i 2 , SIdm'i 2 ,^ t, n s , log A s , A sz , Ej}; and Model 
C, which does not consider any interaction, a seven- 
parameter space {ftbh 2 , 57dm^ 2 , Q, t, n s , log A s , A sz }, 
corresponding to the standard ACDM model. 

The summary of the posterior one-dimensional 
marginalized probabilities is outlined in Fig. 4 and Table 
UJ In Fig. 5 we show the contour confidence intervals for 



4 Special care must be taken because numerical codes, as CAMB, 
use synchronous gauge and fix the residual gauge freedom by 
taking 0dm = 0. In cases in which there are sources in the 9 
evolution equation, like ours, this is not possible to do. 



the marginalized Ej — Ejr space at 0.68 and 0.95 confi- 
dence levels (c.l.). There, the high degeneracy between 
both parameters is shown: while E/i takes values closer 
to zero, E/ also does. It is interesting that nonzero values 
of the interactions (when introduced) are consistent and 
preferred by the considered data at 0.95 c.l. 

We note that the addition of the Ej and E/j inter- 
actions to the theory produce remarkable differences be- 
tween the parameters estimations of Model A and Model 
C. This is more evident for the baryonic matter energy 
density f^/i 2 , as can be observed from Figs. 4, 6, and 7, 
or read directly from Table I. For this reason, we study 
Model B, which only considers the Ej interaction. In 
this case, the tensions between the parameters estima- 
tions are alleviated. The discrepancies are also evident 
in the primordial power spectrum parameters A s and n s , 
which also can be seen from Figs. 4 and 8. Nonetheless, 
these are not alleviated when one considers Model B. 

Instead of using the proton mass as the scale in the 
interactions, we can use an arbitrary associated mass to 
the dark matter or dark fluid particles, m<j. We obtain 
the following constraints at 0.68 c.l. on the ratio E/md 
(we use c = 3 x 10 10 cm/s): 

For the case in which we consider both interactions 
(Model A) 

3.58 x 10- 22 < < 8.68 x 10~ 22 ^/^, (79) 
rrid GeV/c^ 

and 

- 1.94 x 10- 13 < ^££ < -1.05 x 10- 13 7^^. (80) 
m d GeV/c 2 

While for Model B, 

0.22 x 10- 22 < ^ < 1.66 x 10- 22 T^ryyV (81) 
md GeV/c^ 

We only obtained bounds of nid and E in the combina- 
tion E/m^, although it is possible to find constrictions of 
them separately, even to the cross section and the parti- 
cles' thermal velocity, leading to constraints in the a— rrid 
plane. Nevertheless, to do this we have to allow c 2 y^ 0, 
and moreover, make further assumptions on the thermo- 
dynamics of the dark matter (or the dark fluid) [52| • 

Note that the effective thermalized cross section for in- 
teraction iT is a 2 E//, and this is equal to Ej at about a 
redshift z ~ 10 5 . Before this epoch interaction / domi- 
nates, and after that, interaction iT starts to do it. Just 
before recombination, at z ~ 1100, interactions I and 
II are smaller in strength than the Thomson interac- 
tion by about 9 and 5 orders of magnitude, respectively. 
After this time, the ionization fraction x e falls expo- 
nentially and Thomson scattering becomes subdominant 
very quickly. 

In Fig. 9 we plot the CMB power spectrum for the 
mean values estimated by the MCMC analysis. The re- 
ported difference is A = 1(1 + 1) (c, Modcl c - Q). We 
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note that the three models are well inside the error bars 
of the binned measurements from WMAP seven-year 
observations results. Nonetheless, at higher multipoles 
(/ > 1000) the three models start to have notable dis- 
crepancies. The higher differences appear in the region 
1000 < I < 2000, just inside the window in which the pri- 
mordial power spectrum parameters are expected to be 
estimated with higher precision (1000 < I < 3000). For 
I > 3000 secondary anisotropies, mainly the Sunyaev- 
Zeldovich effect, are expected to dominate the CMB spec- 
trum, while for multipoles I < 1000 the spectrum is more 
sensible to the other parameters. These high-/ power 
spectrum multipoles will be probed by the PLANCK mis- 
sion [sH , so we expect to obtain tighter constraints in the 
near future. 

In this work we have not considered the effect that 
interactions I and II could have on big bang nucleosyn- 
thesis. This is because our phenomenological model only 
includes the thermalized cross sections £/ and £// of 
elastic collisions, whose strengths are at least 9 orders of 
magnitude weaker than the Thomson scattering at this 
epoch, and more important, do not annihilate baryons 
and maintain the baryon-to-photon ratio unaltered. Ac- 
cordingly, we expect the effect over this process to be 
quite weak. 



VII. CONCLUSIONS 

In this paper we worked out some properties of the 
so-called dark degeneracy 2l| : the fact that we can only 
measure the total energy momentum tensor of the dark 
sector and any split into different pieces (as in dark mat- 
ter and dark energy) is artificial. Although it could be 
mathematically convenient. 

We start by defining the dark fluid as a barotropic 
fluid with speed of sound equal to zero, as in [29|, l30j 
and, more recently, in (3lj (for similar approaches see 
[H|,|33j]). This is motivated by astrophysical constraints 
in the dark matter besides its cosmological properties. 
Making the speed of sound equal to zero allows dark fluid 
energy density perturbations to grow at all length scales, 
as a cold dark matter component does, while it leaves 
room for a constant pressure different from zero. Astro- 
physical scenarios forbid this pressure to be very high, 
but it well could take values of the cosmological critical 
density today. From the positivity of the energy density 
of the dark fluid it follows that this pressure has to be 
negative. Thus, we conclude that the dark fluid could 
act as dark energy. Then we study the background cos- 
mology for the dark fluid and show that, not surprisingly, 
the same evolution equations as in the ACDM model are 
obtained, leading to a degeneracy between both models. 
We conclude that any collection of fluids for which its to- 
tal equation of state parameter is equal to Eq. (fTg|) will 
lead us to the same result at the background cosmological 
level. 

This is not necessarily the case when we consider linear 



perturbations. In order to preserve the dark degeneracy 
we demand that the dark fluid be a perfect fluid. For a 
barotropic fluid, this last condition implies that the first 
order space-space part of the energy momentum tensor 
is equal to zero. Furthermore, one has to add that the 
energy momentum flux density is equal to that of the 
ACDM. Thus, we have demanded that the energy mo- 
mentum tensor of the dark fluid be equal to the ACDM 
one at first-order in perturbation theory, otherwise the 
degeneracy is broken. 

When considering a collection of multiple interacting 
fluids and demanding that the adiabatic speed of sound 
of the total composed fluid be equal to zero, we obtain 
exactly the perturbation equations of the dark fluid. This 
shows that the degeneracy is present also for more com- 
plicated dark sector schemes. 

From the first five sections, we conclude that there 
exist an infinity of cosmological models that give exactly 
the same observational signatures. Accordingly, it is fun- 
damentally impossible to elucidate the actual structure 
of the dark sector. Thus, it is a matter of taste to prefer 
the ACDM model over any one of the exposed in this 
paper. In fact, to economize it is better to take the one 
fluid choice (the dark fluid) as the correct model. 

The only hope we have to understand the actual na- 
ture of the dark sector, is that it interacts in some way 
with the particles of the standard model. Of course inter- 
actions would in general break the degeneracy, but this is 
not necessarily true. This subject is studied in Sec. VI, 
where we allow the degeneracy of the dark fluid and the 
ACDM model in the homogeneous and isotropic cosmol- 
ogy, and then, we try to break it by adding interactions to 
baryons at first order in perturbations theory. We show 
that a general class of interactions defined by Eqs. (|63|) 
and (|64l) also can be understood as between dark mat- 
ter and baryons in the ACDM model context, leading 
to a new degeneracy, although not with the concordance 
model, but with a ACDM-plus-interactions model. As a 
consequence, these interactions does not help us to favor 
any decomposition of the dark fluid. 

For the latter investigation we have chosen two inde- 
pendent interactions: one that resembles the Thomson 
scattering between photons and baryons, regulated by a 
parameter £/, and a second, density-dependent interac- 
tion, inspired by chameleon theories and parametrized 
by For the analysis we have used a combination of 
the seven-year results of WMAP, the Union 2 data set 
compilation of supernovae measurements, and a prior in 
the Hubble constant of the Hubble Space Telescope. 

When we consider both interactions together, we ob- 
tain that other parameters of the theory differ too much 
from its standard values, up to almost 10% in the case 
of the baryonic density parameter f2&/i 2 . Then, we study 
the case in which £// is equal to zero, allowing S/ to 
vary (Model B). Here, we show that the tension between 
the parameters estimation is reduced. 

The introduction of the two parameterized interactions 
is allowed by current data, and if done, it is remarkable 
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that nonzero values of them are preferred at 0.95 confi- 
dence level. 
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